Preliminary Results with Lattice Covariant Gauge 



ON 



L. Giusti«, M. L. Paciello( 2 \ S. Petrarca^ 2 ' 3 ) , B. Taglienti^ 

1 Boston University - Department of Physics, 590 Commonwealth Avenue 

Boston MA 02215 USA. 
2 INFN, Sezione di Roma 1, P.le A. Moro 2, 1-00185 Roma, Italy. 
3 Dipartimento di Fisica, Universita di Roma " La Sapienza" , 
P.le A. Moro 2, 1-00185 Roma, Italy. 

O 

Q 

psj ! In this poster we present a few preliminary results obtained using our method 

to fix generic covariant gauges on the lattice. We have computed the gluon 
^ • propagator and we have found a sensitive dependence on the gauge parameter. 

As proposed in |], |f , ^ e generic covariant gauge can be defined by replacing 
q . the Landau gauge condition with the following form: 

(N 
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d^A (x) = A(x) + (periodic boundary conditions) (1) 



where A(x) belongs to the Lie algebra of the group. The functional proposed in 
[U in order to fix non-perturbatively the condition (|l|) is 



H[G] = \\d,A G - A|| 2 = J d*xTr [(d,A G - k){d v A° - A)] . (2) 
In fact it has the property 

' D v d v {d,A G - A)l (3) 
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with G(x) = e l9a ^ Ta , where T a are the Gell-Mann matrices. Eq. (|^) shows that 
H[G] is stationary when the eq. (||) is satisfied. Spurious solutions, that can be 
generated by the minimization, do not seem to influence our numerical results. 
On the lattice, in the generic covariant gauge eq. (0), the expectation value of a 
gauge-dependent operator O is obtained by 

(<D) = -J dAe-^I d4xTr ^ J dUO(U G ^)e~^ (4) 

where G(a) is the gauge transformation that minimizes the discretized version of 
the functional (|2|) and the As follow a gaussian distribution of width a. In fact, 
on the lattice, the correct adjustment to the measure is built into the simulation 
recipe and there is no need to compute the Faddeev-Popov determinant. Hence 
the numerical procedure implied by eq. (HI) can be described as follows: 
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• A gauge configuration {U} with periodic boundary conditions according to 
the gauge invariant weight e~ Sw ^ is generated; 

• For each {U} configuration random matrices A(x) belonging to the group 
algebra are extracted according to the gaussian weight of eq. ([§); 

• Given A(x), a numerical algorithm minimizes a discretization of the func- 
tional H[G\. That defines the lattice gauge fixing condition ; 

• The expectation value of the lattice gauge dependent operator is then given 
by the average over the configurations: 

{0) Latt = 1 0{U G ) . (5) 

{conf} 

This is the procedure we will use to compute gauge dependent correlation func- 
tions in a generic covariant gauge. 

Of course, in order to fix the gauge non-perturbatively on the lattice, it is 
necessary to discretize the gauge fixing functional relevant for the gauge condition 
required. The freedom in the choice of the lattice definition of A^, as discussed in 
ref. jl, can be used to build discretizations of the minimizing functional which 
lead to an efficient gauge fixing algorithm. This was the case for the standard 
algorithm of the Landau lattice gauge fixing. It is possible to take advantage of 
the freedom to choose the discretization of the gluon field to find a discretization 
of H[G] ("driven discretization") such that it takes only a local linear dependence 
upon G(x). This aim can be reached by choosing the discretization of each H term 
in order to guarantee the local linear dependence on G(x) instead of following 
a particular definition of A^. Using this idea H[G] can be discretized in the 
following compact form 

Hl(G) = y^TrJ2J G (x)J^(x) (6) 

where 

N(x) = -8/ + E [Ul(x -v) + U u (x)] ; J(x) = N(x) - iagA(x) . (7) 

It is easy to see that locally the functional transforms linearly in G(x) and in 
the continuum limit it goes to the functional (0). The functional (^) is semidefi- 
nite positive and, unlike the Landau case, it is not invariant under global gauge 
transformations. The functional Hl[G] can be minimized using the same iterative 
algorithm adopted in the Landau gauge fixing. In order to follow the convergence 
of the algorithm, two quantities can be monitored as a function of the number of 
iteration steps: the functional Hl[G] itself and 



°h = yY, Tr l A H A Hl where A H (x) = [X H (x) - X%{x) 



Traceless 
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Figure 1: Correlations (AiAi)(t) for different values of a. 



and 



X H (x) 



{^2 U^(x)J(x + fi) + Ul(x - {i)J{x - fj)^j - 8J{x) - 721 



+ 8 (K(x) +K\x)) + N(x)K\x) 



(8) 



Ah is invariant under the transformations A(x) — > A(x) +c like in the continuum. 
9h decreases (not strictly monotonically) reaching zero when Hi\G\ gets constant 
and it signals that the algorithm has converged. In this feasibility study we have 
generated 50, SU(3) thermalized configurations using the Wilson action with 
periodic boundary conditions at (3 = 6.0 for a volume of 8 4 . Then for each value 
of a we have extracted a set of A(x) with a gaussian distribution of width a. We 
have monitored H L [G] and 6 H step by step and the chosen quality is 9 H < 10~ 10 . 
Once the configuration have been gauge rotated we have computed the following 
two point correlation function 



(M)(i) 



w 2 



££Tr(A(x,t)A(y,0)>. 



(9) 



i=l,3 x,y 



In eq. (Q) the standard definition is used for the gluon field. Our result for the 
correlator (AiAi)(t) (H), relevant for the investigation of the QCD gluon sector, 
are shown in Fig. [1|. The relative statistical errors are comparable with the 
Landau gauge case with the same number of configurations. Even with a small 
volume and a small number of configurations, the gauge dependence of the gluon 
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propagator is clearly shown. The a dependence of the gluon propagator shown 
in Fig. |l] does not seem to be re-absorbed by an overall scaling factor. 

This plot shows the feasibility of our procedure to study the gauge dependence 
of physically interesting correlators. 
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